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Abstract 

We calculate the logarithmic moment of the phonon frequency 
spectrum at a single density for 29 monatomic liquids using two meth- 
ods, both suggested by Wallace's theory of liquid dynamics: The first 
method relies on liquid entropy data, the second on neutron scatter- 
ing data in the crystal phase. This theory predicts that for a class of 
elements called "normal melters," including all 29 of these materials, 
the two estimates should closely match, and we find that they agree 
to within a few percent. We also perform the same calculations for 
four "anomalous melters," for which we expect the two estimates to 
differ markedly; we find that they disagree by factors almost up to 
three. From our results we conclude that the liquid entropy estimates 
of the logarithmic moment, applicable both to normal and anomalous 
melters, are trustworthy to a few percent, which makes them reliable 
for use in estimates of various liquid transport coefficients. 



1 Introduction 

Moments of the phonon frequency distribution of systems in the liquid phase 
appear to set the timescales on which certain correlation functions decay, thus 
affecting the values of the transport coefficients to which they correspond by 
Green-Kubo relations [3 EI; so determining these coefficients requires the 
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reliable estimation of the relevant moments. The comprehensive theory of 
monatomic liquids proposed by Wallace (see [S] for a review) allows one to 
estimate one of these moments, the logarithmic moment Bp, from experi- 
mental measurements of the liquid entropy, but for these applications one 
would like greater confidence that this estimate is a reliable one. In the first 
paper on this theory it was shown that the resulting expression for en- 
tropy together with the estimate 0q ~ 0q, where 0q is the corresponding 
moment in the crystal phase (this estimate is explained below) reproduced 
the temperature dependence of the entropy for six elements to an accuracy 
of roughly 2% from melt to three times the melting temperature, supporting 
both the general form for the entropy and the approximation for 0q. How- 
ever, for most materials entropy data well above melting are unavailable, 
and one must ask a different question: If one computes Oq from only a single 
entropy data point near melting, how reliable is the result? We will argue 
that this estimate is very good by considering data for 29 "normal melters" 
(elements like the six in 4j, for which the theory claims Bg ~ Gq is a good 
approximation), showing that 0q computed from entropy data at the den- 
sity of the liquid at melt at 1 bar closely matches 0q computed at the same 
density from neutron scattering data and other sources. The internal consis- 
tency of these two ways of estimating Bq for normal melters, both provided 
by the same theory but utilizing independent sources of data, increases our 
confidence that the entropy estimates of Gq will be trustworthy even for el- 
ements that are not normal melters, and thus aren't expected to satisfy the 
©0 ~ ®o approximation. 

In Section |21 we use this theory to deduce the two ways of estimating 
00 , one of which applies generally and one of which applies only to normal 
melters, and in Section |21 we use experimental data to compute the two 
estimates, finding that they match to 7%. Finally, in Section 0] we consider 
the significance of these results. 

2 Two estimates of 0q for normal melters 

The theory of liquid dynamics that we use here was developed in response 
to two trends in the experimental data for monatomic liquids at melt [3]: 
(1) The nuclear contribution to the specific heat tends to lie near 3k per 
atom in both crystal and liquid phases; and (2) the entropy of melting data 
at constant density naturally divide the liquids into two classes: the "nor- 
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mal melters," for which the entropy clusters around 0.8k per atom, and the 
"anomalous melters," for which the entropy of melting is much higher and 
which undergo significant electronic structure change (e.g., semimetal crys- 
tal to metal liquid) upon melting. The specific heat results suggest that the 
liquid undergoes collective harmonic motion much as the crystal does, and 
the entropy of melting data are naturally interpreted to show that the liquid, 
as opposed to the crystal, moves among w'^ potential valleys, where N is the 
number of atoms and \nw ^ 0.8, thus contributing an extra configurational 
entropy on melting, with the excess for the anomalous melters being due to 
the electronic structure change. These observations led to two hypotheses: 

1. The potential landscape in which the atoms move in the liquid phase is 
dominated by approximately w'^ intersecting nearly-harmonic valleys, 
and the atomic motion is a combination of oscillations in one of these 
valleys and nearly instantaneous transitions from one valley to another 
called transits. Transits are responsible for diffusion. 

2. The valleys fall into three categories: the few crystalline valleys in 
which the system exists as a single crystal, the more numerous sym- 
metric valleys corresponding to configurations that retain some of the 
crystal symmetry, and the overwhelmingly numerically dominant ran- 
dom valleys, which retain no remnant symmetry and which all have the 
same depth and shape (and thus the same phonon frequency spectrum) 
in the large N limit. 

These hypotheses have been extensively tested (as summarized in |H]) and 
can be used to calculate thermodynamic quantities such as entropy per atom; 
in the classical limit, where most liquids reside, the result is (Section 4 of P]) 

S' = k\nw + 3k [ln(T/e[,) + l] + + S^^ + S[, + 5^,,, (1) 

where 5'^„/j, 5*^^^, 5"';, and Sl^^g are the contributions to the entropy from 
anharmonicity, the presence of boundaries of the potential valleys, the ther- 
mal excitation of the electrons, and magnetic effects. (The anharmonic and 
boundary contributions are included in a single term, and the magnetic con- 
tribution is not included at all, in Section 4 of 0, but Smag is discussed on 
pp. 202-203 of 0.) Thus, if one can compute S^^^j, S^^y, S^i, and S^^^g at a 
given density (or argue that they are negligibly small), one can calculate Gq 
from liquid entropy data at that density at a single temperature (or a range 
of temperatures, if data are available). 
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While the theory argues that this means of calculating 0q is valid for all 
monatomic liquids, it also makes another prediction for the value of Gq which 
applies only to the normal melters, namely that 0q is nearly equal to Bq at 
the same density. The prediction comes about as follows (see Subsection 4.4 
of jS]): The entropy per atom of the crystal phase in the classical limit is 
calculated from lattice dynamics to be 

= 3k [ln(r/0S) + 1] + S:^, + SI, + Sl,,^, (2) 

where 5"^^;^, 5"^;, and S^^g represent the contributions from anharmonicity, 
the electrons, and magnetic effects, respectively, so the entropy of melting at 
constant density AS' is given by 

AS = S'(TJ-S=(TJ 

Sl^y is omitted here because boundary effects are non-negligible only at tem- 
peratures well above melting. For a normal melter, the lack of significant 
change in electronic structure suggests that S'^ ~ S^i, and assuming that an- 
harmonic and magnetic effects are small (or at least comparable in the two 
phases), we predict AS ~ /c In -(- 3/cln(0Q/0Q). Experiment shows that for 
the normal melters AS = 0.8A; with a small scatter, which strongly suggests 
that Inw = 0.8 (since w is the only available parameter that is not strongly 
material-dependent) and that Gq ~ Gg, with anharmonicity and small dif- 
ferences between the two values of Gq accounting for the scatter. Thus one 
should be able to estimate Gq for normal melters reasonably closely by de- 
termining Gq from appropriate data. 

If this interpretation of the entropy of melting data is correct, these two 
estimates of Gq for normal melters should approximately match, which will 
in turn increase our confidence in estimates of Gq from liquid entropy data 
for any liquid, normal or anomalous. Determining whether they do in fact 
match is the subject of the next Section. 

3 Comparison of experimental estimates 

Let pim be the density of the liquid at melt at 1 bar (with the exception of 
Ar; see the Appendix). We will compare GQ(p/m) and Go(/9/m) for 29 normal 
melters for which accurate data are available. 
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To compute QQ^pim), we consult neutron scattering or other data which 
give us 0o(Pmeas), whcrc Pmeas IS souie density not too far from pi^, and 7*, 
the high-temperature Griineisen parameter. From these two quantities one 
can calculate Oo(pzm) from the relation 

Ql{Plm) = e'oiPmeas)(-^y , (4) 

\Pmeas J 

which is highly accurate as long as \pim — pmeas\ is not too large. 

To compute Oo(pim) from liquid entropy data, we first define the harmonic 
contribution to the entropy S\^^^^ by 



Sharm = klliW + 'ik 




(5) 



where is the quadratic moment of the liquid phonon distribution. This is 
the same as the first two terms in Eq. for the classical limit of the liquid 
entropy, carried to the next higher order term in T~^. Then if Sl^^p^ is the 
experimental value of the liquid entropy at the melting temperature Tm at 1 
bar, 

^expt ~ ^harmiPlm, ^m) + >S'^n?t ~^ ^el ~^ ^mag (6) 

(again, Sl^y is negligible at melt), so to determine 5'^^^™ (and thus Gg) we 
need 63 and the remaining three terms on the right hand side of Eq. (jHl). 
The standard approximation B2 = c^^^Oq is perfectly reasonable here, since 
the 02 term contributes less than 1% to the entropy at the values of Tm 
we consider, and we argue that the anharmonic and magnetic terms in the 
liquid at melt will be roughly equal to their counterparts in the crystal just 
before melt. (Recall that these two terms can be calculated for the crystal, 
because we have independent data for Gq and thus we can determine the 
harmonic contribution to the crystal entropy directly; once the electronic 
term is also found from theory, the anharmonic and magnetic parts are simply 
what's left.) The electronic term can be calculated different ways for different 
materials; see the Appendix for a discussion of how 5'^ was calculated for 
the materials here. 

These calculations were performed for 29 normal melters, using the data 
collected in the Appendix, and the results are shown in Table ^ In Figure 
HI 0o(Pim) is plotted against Qo{pim) for all 29 elements, along with the hne 
0o(p«m) = 0o(p«m) for comparison. First, we note that all 29 points cluster 
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Element 9^ (K) 0[, (K) &y&'o 



Cs 


26.4 


27.2 


1.03 


Cd 


86.6 


84 


0.97 


Tl 


(54.5) 


55.2 


(1.01) 


Li 


252 


251 


1.00 


Na 


102.2 


100 


0.98 


K 


61.9 


59.7 


0.96 


Rb 


37.7 


37.3 


0.99 


Mg 


200.4 


181 


0.91 


Sr 


(76.3) 


70.7 


(0.93) 


Cu 


176.9 


172 


0.97 


Ag 


113.6 


118 


1.04 


Au 


91.5 


97.4 


1.06 


Zn 


129.3 


124 


0.96 


Hg 


56.4 


52.8 


0.94 


Al 


209.0 


197 


0.94 


In 


74.4 


75.0 


1.01 


Pb 


51.4 


54.1 


1.05 


Ar 


59.1 


67.8 


1.15 


V 


212 


191 


0.90 


Nb 


169 


161 


0.95 


Ta 


134 


121 


0.90 


Pt 


119 


123 


1.03 


Pd 


145 


(146) 


(1.01) 


Ni 


218 


209 


0.96 


Cr 


(276) 


(283) 


(1.03) 


Mo 


235 


(201) 


(0.86) 


W 


188 


(165) 


(0.88) 


Ti 


(197) 


204 


(1.04) 


Zr 


(141) 


127 


(0.90) 



Table 1: 0q(p;^) and QQ{pim) ioi 29 normal melters, together with their 
ratios. Less reliable numbers are shown in parentheses; see the Appendix for 
a discussion of the rehabihty of data. 
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Figure 1: Qo(pim) versus 0o(p;m) for the 29 elements in Table [H The seven 
elements indicated by solid circles, all but one of which are transition metals, 
deviate systematically from the line QQ{pim) = ^oiPim)] we believe this is due 
either to these materials being weakly anomalous melters or to uncertainties 
in the entropy data. 



around the line to good approximation; second, we note that most of the 
scatter is due to the systematic deviation of the seven materials indicated by 
solid circles, all but one of which (Mg) are transition metals. This could be 
due either to the larger uncertainty in the experimental entropy results for 
the transition metals or to the fact that some of these materials are weakly 
anomalous melters (the entropy of melting data for Mo and W particularly 
suggest this |5|). Our data for all 29 materials show 

-f = 0.978 ± 0.063 (7) 

at pim, or equivalently that the rms deviation of Qq/Qq from 1 is 0.067 at 
pim- As we explain in the Appendix, given the quality of our data we expect 
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Element 6^ (K) e[, (K) elj/Gfj 



Ga 


169 


100 


0.59 


Si 


436 


162 


0.37 


Ge 


255 


91 


0.36 


Sn 


94.3 


74 


0.78 



Table 2: QQ{pim) and 0o(pim) for four anomalous melters, together with their 
ratios. As expected, none of the ratios are near unity, reflecting the changes 
in the interatomic forces (and thus phonon frequencies) that occur at melt 
for these materials. 

roughly 3% (at most 4%) errors in this ratio; so as the theory predicts, 
the two are equal to a high degree of approximation, and about half of the 
difference is accounted for by experimental error. 

In contrast, if we consider analogous calculations for four anomalous 
melters, we find the results shown in Table |21 The change in electronic 
structure at melt results in great changes in the interatomic forces, resulting 
in large changes in the phonon frequencies. For these materials, the approx- 
imation 9q = 9q would be grossly inaccurate; but as these elements are 
nearly-free-electron metals in the liquid phase, we expect the calculation of 
0Q from entropy data to be quite reliable. 

4 Conclusions 

The theory used here gives us two ways to calculate 6q, one using liquid 
entropy data and valid for all monatomic liquids, and the other using crystal 
neutron scattering data and applicable only to "normal melters," for which 
the crystal and liquid values of Go as a function of density are expected to 
agree approximately. This expectation is amply fulfilled by data from 29 
normal melters, which show Bq and Gq agreeing within a few percent. (Our 
expectation that 0q and 9q of anomalous melters should strongly disagree 
is also verified.) This work nicely complements calculations by Wallace and 
Clements for pseudopotential sodium, in which they found at p = 0.925 
g/cm^ 

e^^'' = 99.65 K and = 98.7 K, (8) 

which agrees very well with the experimental results for Na from Table d 
These results, supporting the internal consistency of this theory, suggest that 
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estimates of 0q from liquid entropy data for anomalous melters, such as those 
in Table |2l will be equally trustworthy; coupled with the theoretical result 
given in Eq. (jH)), they also show that theoretical methods have advanced to 
the point that a priori calculations of 0q, which will be important for the 
computation of transport coefficients in the liquid, are also reliable. 

Appendix: Data and Sources 

All of the data used to calculate Oo(pim) and Oo(pzm) for the 33 elements 
we consider are recorded in Table 121 The entry in the "Structure" column 
indicates the structure of the crystal at melt; Pmeas is the density at which 
6q was measured; 7* is the high-T Griineisen parameter; pi^ is the density 
of the liquid at melt at 1 bar (except for Ar; see below); is the corre- 
sponding melting temperature; S^^^^ is the experimentally measured entropy 
of the liquid at melt; S'g^, and are the additional contributions 

to the total entropy from Eq. ((H) {Sl^y is negligible at T^); and S'^^rm is 
the harmonic contribution to the entropy determined from the other entropy 
data using Eq. All densities are in g/cm^, all temperatures are in K, and 
all entropies are in fc/atom. Data in parentheses are generally less reliable. 

Most of these data are taken from information collected in various tables 
in jH]. For example, the structure information is taken from Table 19.1, pp. 
191ff. Values of 6o(pmeas) for materials in the first group were determined 
using crystal entropy data collected in Problem 19.8, p. 207 (neutron scat- 
tering data is either unavailable or was taken in a phase other than the one 
from which the crystal melts); the second, third, and fourth groups were de- 
termined using neutron scattering data recorded in Table 15.1, p. 153; and 
the final two elements were handled by taking high-T neutron scattering data 
from and inferring Bg at (see Figure 9 of j?!). These renormalized Gq 
include anharmonic contributions, so S^^f^ was neglected for these elements. 
(The fact that Gg and Gq are so close for these two materials [see Table [T] 
suggests that the crystal and liquid phases have approximately the same an- 
harmonicity; this is also suggested by the specific heat data from Figure 1 of 
^.) The corresponding values of pmeas for all materials were calculated from 
the experimental temperatures using extensive tables of density-temperature 
data at 1 bar for the elements. 

The 7* data are also found in Table 19.1 of |3] (except for the few val- 
ues in square brackets, which are purely empirical estimates where data are 
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Element 


Strueture 


Pmeas 


Oq {.Pmeas ) 


7* 


Plm 


T 


0/ 

expt 


0/ 


0/ 

anh 


ol 
maq 


ol 
harm 


Cs 


bee 


1.91 


27.5 


1.14 


1.84 


301.6 


11.11 


0.09 








11.02 




hep 


o.OO 


iUo 


9 Q 
z.o 


S 09 
O.UZ 


t^QA 9 

oy4.z 


Q 7f; 
y. ( 


n n7 

U.U ( 


n 
u 


n 
u 


fiO 

y .oy 


1 1 


bee 


ii.O 


Oo.o 


f9l 


1 1 99 
ii. zZ 


c^77 


1 n OQ 
iu.yo 


n nn 
u.uy 


n 
u 


n 
u 


1 fi/i 

1U.04 




ortho 


c; Ol 
o.y i 


1 fi9 


1 c; 

i.O 


R no 
o.uy 


Qn9 
oUz.y 


7 1 s 
( .io 


n n/1 
U.U4 


n 
u 


u 


7 1/1 
( .i4 


ol 


clia 


Z.o4 


/I 91 
4Zi 


(0.5j 


9 c;i 
z.Oi 


1 Rfi7 
iOo ( 


1 1 nc; 
ii.UO 


n 99 
U.zz 


n 
U 


n 
U 


1 fiQ 
iU. 00 




dia 


o.o/ 


9/1 c; 
Z40 


U.o 


c; fin 
O.OU 


1911 

izii 


1 1 79 

ii. ( z 


n 1 7 
U. i ( 


n 
U 


n 
U 


1 1 

ii.OO 


on 


DCt 


( .OU 


1 DQ /I 


9 9 
z.z 


7 nn 
( .UU 


KHK 1 
OUO. i 


y.oo 


n no 
U.uy 


n 
U 


n 
U 


c;7 
y.o ( 


Lil 


bee 


U.04u 


ZuO.O 


U.oo 


n c;i 

U.OiO 


/I 7 

40o. ( 


0.00 


n n/1 
U.U4 


n 
u 


n 
u 


R9 
o.oz 


Mo 


bee 


i.UUO 




1 9/1 
i.z4 


u.yzo 


Q71 n 
( i.U 


7 7fi 
( . ( 


n nc; 

U.UO 


n 
U 


n 
U 


7 7Q 
1.(0 


IV 


bee 


n on/i 
U.yU4 


Do.y 


1 9/1 
i.z4 


u.ozy 


QQR /I 


nfi 
y.uo 


n n7 
U.U ( 


n 
U 


n 
U 


s 00 
o.yy 


rvD 


bee 


i.DiO 


/I9 9 
4Z.Z 


1 9fi 


1 /1 7Q 

i.4 ( y 


'ii 9 R 
OiZ.O 


1 n 9R 

iU.ZO 


n ns 

U.Uo 


n 
u 


n 
u 


10 18 
iU. io 


A/Trr 

Mg 


hep 


1 7/1 


990 A 

zzy.4 


1.0 


i.oy 


099 

yzz 


S fil 

o.oi 


n 1 9 
u. iz 


n 
U 


n 
U 


o.oy 


or 


bee 


9 R 
Z.O 


oU.U 




9 /I s 
Z.4o 


1 n/1 9 

iU4Z 


1 9 1 n 

iZ. iU 


n 9^? 
u.zo 


n 
u 


n 
u 


1 1 87 


nn 

uu 


fee 


y.uio 


99c; Q 


9 n9 

Z.UZ 


o.UU 


ioOo 


1 n no 
iU.uy 


n no 
u.uy 


n 
u 


n 
u 


1 00 
iU.UU 


Ag 


fee 


1 n /lo 
iu.4y 


1 e;n 1 
iou.i 


9/19 
Z.4Z 


y.oO 


1 9Qc; 
izOO 


1 n OA 
iU.y4 


010 

U. iU 


n 
U 


n 
U 


1 S/l 
iU.o4 


An 

All 


fee 


1 O 97 

iy./ ( 


1 9/1 e; 


9 Qc; 
z.yo 


1 7 Qfi 
i ( .OD 


1 QQfi 

iOoo 


1 1 77 

ii. M 


oil 

U. i i 


n 
U 


n 
U 


1 1 fiR 
ii.OO 


Zn 


hep 


7.270 


161 


2.2 


6.58 


692.7 


9.04 


0.07 








8.97 


Hg 


rhomb 


14.46 


64.7 


2.5 


13.69 


234.3 


8.31 


0.03 








8.28 


Al 


fee 


2.731 


283.5 


2.25 


2.385 


933.5 


8.58 


0.11 








8.47 


In 


bet 


7.43 


85.3 


2.4 


7.02 


429.8 


9.11 


0.07 








9.04 


PK 


fee 


A.A..OO 


f\A 1 

U^. -L 


9 fi^ 
z.uo 


in f,"^ 

lu.uo 


ouu.o 


1 1 1 'i 
J- J- . 10 


n 1 1 


n 
u 


n 
u 


1 1 09 
J- J- .uz 


At- 

Ar 


fee 


1 77nc; 
i. ( ( Uo 


c;n 1 
oy.i 


9 R 
z.o 


1 77nc; 
i. ( ( UO 


91 

ziy 


7 QQ 
( .00 


n 
U 


n 
U 


n 
U 


7 QQ 
( .00 


V 


bee 


0. iZ 


zou 


1 

i.O 


o.4y 


99n9 

zzuz 


(^o 1 9^ 

iZ. iZ J 


08 
u.yo 


n 

U 


n 

U 


(^ii.i4j 


IN D 


bee 


o.OO 


1 Ofi 

iyo 


1 

i.O 


7 7fi 


97/1 A 
Z ( 44 


1 Q 9/1 

io./4 


O/l 

u.y4 


n 
u 


n 
u 


1 9 QO 

iz.oU 


Id 


bee 


lu.uo 


1 fi9 
1 uz 


l.U 




OZb/O 


1 A 1^ 


1 (XX 

l.UO 


n 
u 


n 
u 


1 79 
10. 1 z 


JT I 


fee 


91 ^f, 


LkJO.O 


9 fi 
z.u 


1 Q 1 


zu^o 


1 'i DQ 
10. uy 


n 8^1 

u.oo 


n 
u 


n 
u 


19 94 


Pd 


fee 


12.07 


200 


2.3 


10.5 


1827 


(12.20) 


0.82 








(11.38) 


Ni 


fee 


8.91 


275 


1.9 


7.89 


1728 


11.54 


1.09 





0.31 


10.14 




bee 


7 1 O 

( .iy 


QQfi 
000 


(1.5j 


[v.Ao ) 


91 QQ 
ziOO 


1 1 RO 

ii.oy 


^'O fi/l ^ 
(^U.o4 ) 


[v.vo ) 


n Q/i 
U.o4 


(^y.oo ) 


Mo 


bee 


10.22 


273 


1.6 


9.3 


2896 


(13.26) 


0.67 


(0.79) 





(11.80) 


W 


bee 


19.26 


225 


1.6 


17.2 


3695 


(14.52) 


(0.75) 


(0.65) 





(13.12) 


Ti 


bee 


4.24 


202 


[1] 


4.14 


1945 


11.83 


1.26 


NA 





10.57 


Zr 


bee 


6.25 


147 


[1] 


(6.0) 


2128 


13.39 


1.14 


NA 





12.25 



Table 3: The data used to ealeulate 0o(p/m) and QKpim) for the 33 elements 
in Tables^ and m The units of the data, sourees and meanings of the various 
headings and symbols are diseussed in the text of the Appendix. 
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unavailable), and the pim, Tm, and S^.^^^ data are from Table 21.1, pp. 220ff 
(again, except for Ar; see the next paragraph). The electronic entropy S^i 
was determined for the elements in the first two groups using a free-electron 
model; for Ar is zero because Ar is not a metal; and results for the last 
two groups were computed using density of states calculations by Eriksson 
et al. 0- S^n/i was computed, as mentioned in the main text, by assuming 
it is equal to S^^f^, which is negligible for all but the last five materials in 
the list (see Tables 19.2 and 19.3, pp. 201 and 203 of p^). Anharmonicity in 
the last two is taken into account by the renormalized Qq (see the previous 
paragraph), and for Cr, Mo, and W, S'^„^ (and thus 5'^„^) was estimated by 
combining anharmonic entropy calculations from [7 and Table 19.3 of jH] to 
extrapolate 5"^^;^ to the densities in question. Finally, Smag for ferromagnetic 
Ni and antiferromagnetic Cr were also taken from Table 19.3. 

Special treatment was required for Ar because if pmeas is the density 
at which the neutron scattering data were taken and pim is the density of 
the hquid at melt at 1 bar, \pmeas — Pim\ is far too large for Eq. (jU to be 
reliable. Thus we decided to set pim = Pmeas, as shown in Table El and used 
very extensive tabulated data of the Ar melt curve to find the pressure such 
that Ar would melt at the required density; from this we were able to infer 
the tabulated values for and S^^p^, and from that point the calculations 
proceeded as with the other materials. 

Finally, a few words are in order concerning the quality and accuracy 
of our data. On the crystal side, by comparing phonon moments computed 
from Born- von Karman fits to various different sources of neutron scattering 
data, one sees that determinations of 6q (pmeas) are typically accurate to 
0.5%; given that measurements of 7* typically have 10% error or less, an 
additional error of 2% or so is introduced by the computation of QQ^pim)- 
These errors are independent, so the final error in OQ(p/m) lies between 2 and 
2.5%. For the liquid, comparison of liquid entropy data from separate high- 
quality compendia, such as [IHIIIII among others, indicates errors of ±0.5% 
for most values of S^^p^. Various different electronic structure calculations 
of S'g; vary by 5-10% at most; these errors will affect our entropy results on 
about the 0.1% level for the materials in the first two groups and 1% for 
those in the last two groups. Finally, we estimate the error in S^^f^ + Sl^^g 
to be about 0.5% of S^^p^, so the harmonic part of the liquid entropy S^^j.^^ 
(and thus Gg) will be trustworthy to between 1% and 1.5%. Thus the error 
in the ratio is at most 3 to 4%; this is the source of the error estimate given 
in Section 121 
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Eight materials which have been studied in connection with melting pre- 
viously but are not included in this analysis are Kr and Xe, for which there 
are no entropy data along the melt curve at sufficiently high compression 
(so they cannot be treated as Ar was); and Fe, Th, Sb, Bi, Ba, and Ca, for 
which there are no reliable calculations or estimates of S^i- This is a signifi- 
cant problem with Fe and Th, for which S^i is expected to be large; it is not 
such a difficulty with the remaining four, but when combined with various 
other problems with their data (for example, neutron scattering was done on 
Ca in the fee phase while the crystal actually melts from bcc), it is proper 
to neglect them here as well. 
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